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ABSTRACT 

We study the binary fraction of the globular cluster MIO (NGC 6254) as a 
function of radius from the cluster core to the outskirts, by means of a quan- 
titative analysis of the color distribution of stars relative to the fiducial main 
sequence. By taking advantage of two data-sets, acquired with the Advanced 
Camera for Survey and the Wide Field Planetary Camera 2 on board the Hubble 
Space Telescope, we have studied both the core and the external regions of the 
cluster. The binary fraction is found to decrease from ~ 14% within the core, 
to ~ 1.5% in a region between 1 and 2 half- mass radii from the cluster centre. 
Such a trend and the derived values are in agreement with previous results ob- 
tained in clusters of comparable total magnitude. The estimated binary fraction 
is sufficient to account for the suppression of mass segregation observed in MIO, 
without any need to invoke the presence of an intermediate-mass black hole in 
its centre. 

Subject headings: binaries: general; globular clusters: individual (MIO, NGC6254) 



^ Based on observations collected with the NASA/ESA HST^ obtained at the Space Telescope Science 
Institute, which is operated by AURA, Inc., under NASA contract NAS5-26555. 
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INTRODUCTION 



The binary fraction is an essential component in the formation and evolution of dynam- 
ically active systems, like globular clusters (GCs). In such dense environments, where stellar 
gravitational interactions are very frequent, binaries can exert a significant influence on both 
the dynamical evolution of the system and the properties of its stellar populations. 

Being, on average, more massive than the other stars, binaries tend to sink into the 
highly crowded cluster centers, because of equipartition. The characteristic timescale of this 
process (the relaxation time) depends on the cluster structure and can be even longer than 
a Hubble time in the outskirts. Hence, in the outer regions of GCs we essentially expect to 
observe primordial binary systems, i.e., binaries created as part of the star formation process 
and evolving undisturbed. In the cluster core, on the other hand, a variety of dynamical 
processes (exchange interactions, three-body encounters, tidal captures, etc.) can take place, 
with competin g effects on the binary population: binaries can be destroyed, created or just 



modified (e.g. Hut et al. 1992h. with rela t ive efficiencies that still are a matter of debate 



in th e literature (e.g. Ilvanova et al.l l2005l : iHurley et al.l 120071 : ISoUimal l2008l : iFregeau et al. 



20091 ). In general, however, the gravitational encounters occurring in the core tend to make 
the binaries harder (more tight ly bound) , thus providi ng a central energy source able to slow 
down the c luster core collapse ( iGoodman fc Hutlll989l ). Following the N-body simulations of 
Gill et al.l ( 120081 ) . this energy source could also suppress the mass segregation process, with 
a detectable effect on the radial behavior of the mass function of main-sequence (MS) stars. 
Since the same effect could alternatively be due to a cen tral intermediate- mass black hole 
(IMBH; see also iPasquato et al.ll2009l : iBeccari et al.ll2010l . herefater BIO), if we can measure 
the fraction of binaries, then we can say whether or not we need an IMBH to explain the low 
level of mass segregation that has been observed. Hence, the empirical estimate of the binary 
fraction in a sample of GCs representative of different environments is a prime ingredient for 
dynamical models, which help us understand the internal cluster dynamics. 

The knowledge of the binary fraction is also crucial for understanding the properties of 
puzzling objects like blue stragglers, millisecond pulsars a nd cataclysmic variables, which are 
all thought to be the by-products of bin ary evolution (e.g. jMcCrea lll964j : lRomani et al.lll987l : 
Ferraro et al. Il200ll : iLeigh et al.ll201ll and references therein). In particular, the analysis of 
the bimodal radia l distribution of blue strag glers observed in a number of GCs (e.g., Ferraro 
et al. 1997, 2004; iDalessandro et al. Il2008[ ) suggests that a non- negligible fraction of these 
stars is generated by primordial binaries, which still orbit in isolation in the cluster outskirts 
and produce the observed rising branch of the distribution (Mapelli et al. 2004, 2006; Lanzoni 
et al. 2007a,b). The interpretation of the double blue straggler sequence rece ntly discovered 



in the core of M30 also requires a significant fraction of primordial binaries (IFerraro et al. 
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20091). 



Despite its implications, however, the binary fraction in GCs still remains badly con- 



monly used for its estimate are: radial velocity variability survevs fe.g. iPrvor et al. 



Latham 1996 : Albrow et al. searches for eclipsing binaries (e.g. lMateo lll996[ICote et aL 



(e.g. 


Prvor et al. 


1989 




Mateo 


1996:ICote et al. 



19961 ) , and the stud y of the distribution of stars along the cluster MS in color-magnitude dia- 



grams (CMDs; e.g.. 



Romani et al. 111991 



Bolte 



1992 



a 



Rubenstein fc BailvJl997l: iBellazzini et al 



2002l : IClark. Sandquist fc Bolte 



2004 



Zhao fc Baily: 



2005 



SoUima et al. 



2007 



Milone et al. 



20081 ) . The first two methods rely on the detection of individual binary systems in a given 
range of periods and mass ratios. Hence, the nature of these methods leads to intrinsic 
observational biases and a low detection efficiency. The latter approach relies on the simple 
fact that, since the flux of unresolved binaries is equal to the sum of the fluxes of the two 
components, the binaries composed by MS companions are shifted towards brighter magni- 
tudes with respect to the single-star MS. This technique has the advantage of being more 
efficient and detecting binary systems regardless of their orbital periods and inclinations. 

For the present paper we used this latter technique to estimate the binary fraction in 
the core and the outskirts of MIO (NGC 6254). T his is an "ord inary", dynamically relaxed 

2010 edition) , central mass 



arris 



1996 



GC, with absolute visual magnitude My = —7.48 (Mi 
density logpo = 3.8 (po being in units of Mfn/pc^; 
relaxation time th ~ 0.8 Gyr (Harris 1996; see also iMcLaughlin &: van der Marelll2005l ). The 



Prvor fc MevlanI 119931) . and hal f-mass 



deep and high-quality photometry that BIO obtained for both the centre and beyond the half- 
mass radius, allowed them to study the cluster mass function at different radial distances. 
The resulting mass-segregation profile is moderately flattened and can be explained by the 
presence of either an IMBH of ~ IO^Mq, or a population o f binaries with an initial fraction 
of 3-5% (BlO). Hence, within the framework proposed by iGill et al.l ( 120081 ). any empirical 
constraint on the binary content in this system would allow us to assess the possible presence 
of a central IMBH. In addition, it will provide precious clues and constraints that would 
be useful for a robust interpretation of the properties the blue-straggler-population in this 
cluster (Emanuele Dalessandro et al. 2011 in preparation). 

The paper is organized as follows. The used data-sets are presented in Section 2. The 
method adopted to estimate the binary fraction is outlined in Section 3. The results and the 
discussion are presented in Section 4 and 5. 
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2. THE DATA 

The data-set used in the present work (the same as in BIO) consists of a sample of 
4 X 90sec images acquired in the F606W {V) and 4 x 90sec images in F814W filters obtained 
with the Advanced Camera for Surveys (ACS; GO-10775, PI: Sarajedini), complemented 
with 2 X llOOsec and 2 x 1200sec images in F606W and 2 x llOOsec and 2 x 1200sec images 
obtained with the Field Planetary Camera 2 (WFPC2; GO-6113, PI: Paresce) on board the 
Hubble Space Telescope. The ACS data-set samples the cluster central regions, while the 
WFPC2 one covers an area located between one and two half-mass radii (see Figured]). The 
detailed description of the data reduction, photometric calibratiorj^ and astrometric solution 
procedures is given in BIO and the sample used in the present analysis contains the same 
"bona-fide" stars selected on the basis of the quality of the point-sp read-function fitting (as 



measured by the DAOPHOTII sharpness parameter; IStetsod 119871 ). The CMDs of the two 
data-sets are shown in Figure [2l Stars brighter than / = 16 and / = 19.5 are saturated 
in the ACS and WFPC2 samples, respectively. BIO also performed a detailed photometric 
completeness study, based on artificial star experiments, where artificial stars were added to 
the original FLT frames and the whole data reduction procedure was repeated. From the 
resulting catalogue, listing the input and output positions and magnitudes for more than 
500,000 artificial stars, BIO estimated that the photometric completeness drops below 50% 
at J ~ 22.5 in the innermost region of the cluster, and at / ~ 25 for the WFPC2 data-set. 



3. THE ANALYSIS 



In orde r to estimate th e bina ry fr action of MIO, we fo llowed the method extensively 
described in iBellazzini et al.l (120021 ) and ISollima et al.l (120071 hereafter S07). The basic idea 
is that the magnitude of a binary system corresponds to the luminosity of the primary (more 
massive) star, increased by the contribution of the companion of an amount that depends 
on the mass ratio of the two components [q = M2/M1). In fact, since the stars along the 
MS follow a mass-luminosity relation, the luminosity of a binary can be written in terms of 
the mass ratio of the two components. By definition, < g < 1 and for g = 1 (equal-mass 
binary) the system appears ~ 0.75 magnitudes brighter than the single component, while 
the luminosity enhancement decreases for decreasing q. The spanning of all the possible 
values of q, at different magnitudes of the primary component, produces a broadening of the 



^Thc instrume ntal V and / magnit udes have been calibr a ted t o the VEGAMAG system by following 
the prescription of ISirianni et al.l (|2005l ) and lHoltzman et al. I (|l995l ) for the ACS and the WFPC2 samples, 
respectively. 
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single-star MS, to its bright- and red-hand side. In principle, the ratio between the number 
of stars lying on the red side of the single-star MS and the total number of stars observed 
along the "broadened MS" provides the cluster binary fraction. In practice, depending on 
the photometric error of the data, a minimum value of the mass ratio (qmin) exists below 
which it is impossible to observationally distinguish a binary system from a single MS star. 
Moreover, it is necessary to take into account a number of effects, like stellar blends and the 
contamination by foreground/background field stars, which can add spurious sources in the 
CMD. 

Indeed, chance superpositions of two stars (blends) can produce a luminosity enhance- 
ment that mimics the magnitude shift characteristic of a genuine binary system. In order to 
correct for this effect we analyzed the distribution of the residuals between the input and the 
output magnitudes of the artificial star catalogue built by BIO for the completness study (see 
previous section). From the asymmetry of the distribution (which is skewed toward brighter 
output magnitudes because of the blending between artificial and real stars) we estimated 
that the percentage of blended sources that would mimick binary systems with q > qmin, 
varies from ~ 6% in the core, to less than 0.2% in the external regions. 

BIO also estimated the Galactic field contamination in the direction of MIO, finding 
that it is very low: even in the worst case (the WFPC2 data-set), where the number of 
cluster sources is small, the field stars are just ~ 3% of the total sample. Despite such a low 
value, for a proper measurement of the binary fraction we performed a detailed study of the 



field c ontamination as a function of the magnitude. From the Galaxy model of iRobin et al. 



we retrieved a catalogue covering an area of 0.5 deg^ in the direction of MIO, and 



we randomly extracted two sub-samples of synthetic stars, scaled to the fields of view of 
the ACS and WFPC2 data-sets. Their magnitudes were converted from the Joh n son t o 



the VEGAMAG photometric system adopting the prescriptions of ISirianni et al.l ( l2005l ). 



Finally, by exploiting the artificial-star catalogue used for the completeness studjo (Sect. 
[2]), we obtained a catalogue of synthetic field stars that includes the observational biases 
(incompleteness and blending), for both the ACS and WFPC2 data-sets. 

Once all the contaminant effects are taken into account, the binary fraction was esti- 
mated as the number of stars in the "binary population" divided by the total number of 
stars, i.e., binaries plus genuine, single, MS stars (hereafter the "MS population") . The "MS 



•^publicly available at |http://inodel.obs-besancon.fr/ 

"^For each considered synthetic field star, we randomly extracted an artificial object with similar magnitude 
(A/ < 0.1) and we assigned the shifts between its input and output magnitudes to the field star, in a effort 
to mimick and take into account the effects of completeness and blending. 
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population" is defined as the set of stars liaving a color difference from the MS ridge line 
(MSRL) smaller than three times the typical photometric error at that magnitude level (see 
Figure E]). The operational definition of the "binary population" is given in Sects. 14.11 and 



4. RESULTS 

The high photometric quality and the spatial coverage of the data-sets previously de- 
scribed allowed us to study the binary fraction at different distances from the cluster centre. 
In particular, here we have defined three concentric annuli bounded by the core radius and 
the half-mass radius. The adopted centre of gravity and structural parameters have been 
recently determined from resolved star counts (Dalessandro et al. 2011): the coordinates of 
the centre are aj2ooo = 16*^ 57™ 8.92*^, 5j2ooo = —4° 5' 58.07"; the core, half-mass and tidal 
radii are = 48", = 147", and = 19.3', respectively. This center is located at ~ 3.5" 
North- West from the one quoted by Goldsbury et al 2010, a difference that has no impact 
on the following analysis and the obtained results. Hence, the first two radial bins (r < Tc 
and Tc < r < Vh) are sampled by the ACS data-set, while the third one (r > Vh) is cov- 
ered by the WFPC2 data (see Fig. [1]). Since the two data-sets have different saturation 
and completeness levels (see Sect. [2]), we perfomed the analysis in two different magnitude 
ranges: the adopted cuts along the MSRL are 18.8 < I < 21.5 for the ACS sample, and 
20.3 < J < 23 for the WFPC2 one (see Figs. [2] and [3]). These intervals define what we call 
the magnitude range" of the two data-sets. Then, with the aim of having an interval of 
magnitudes in common between the two samples where to directly compare the computed 
binary fractions, we considered three magnitude sub-ranges defined as follows: a bright 
range" corresponding to 18.8 < I < 20.3, an ^intermediate range" at 20.3 < I < 21.5, and 
a faint range" at 21.5 < / < 23 (all the quoted magnitude values are measured along the 
MSRL). As is apparent from Figs. [2] and [3l the bright range is probed only by the ACS 
data-set, the faint range is found only in the WFPC2 sample, while the intermediate range 
is in common between the two. 



4.1. The minimum binary fraction 

We first estimated the minimum binary fraction {^mm), which is the fraction of binary 
systems with a mass ratio g^m large enough to make them clearly distinguishable from 
the single-star MS. It is clear that the value of g^m depends directly on the photometric 
errors and represents only a sub-sample of the whole population of binaries, but it has 
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the advantage of being a purely observational quantity. In this case we define the "binary 
population" as the set of stars located in the CMD between the following boundaries (see 
gray region in Fig. [3]): the left-hand boundary is the line corresponding to a color difference 
from the MSRL equal to three times the photometric error at any magnitude level (right 
dashed line); the right-hand boundary is the line at a color difference from the equal-mass 
binary sequence equal to three times the photometric error; the upper and lower boundaries 
are set by the largest and the smallest primary mass (corresponding to the quoted bright 
and faint cuts of the various magnitude ranges along the MSRL), combined with all the 
possible mass ratios. In other words, the "binary population" includes all binary systems 
with primary mass set by the considered magnitude ranges and with g^jn < Q' < 1, also 
taking into account the effect of photometric errors. 

For each of the considered radial bins and magnitude ranges we estimated the minimum 
binary fraction by performing all the steps described in Sect. [3] and, in much more detail, 
in S07. The results are presented in Table [H As is apparent, ^min monotonically decreases 
from the center to the outskirts, in agreement with previous findings and with theoretical 
predictions (see Sect. In the full magnitude range, such a radial variation ranges from 
~ 6% at r < Tc, to ~ 1% at r > r/^. There also seems to be a trend with magnitude, 
especially in the central bin, where C,min varies from ~ 8% in the bright range, to ~ 5% in 
the intermediate one. 

However, since the photometric error depends on magnitude, the value of qmin changes 
in the considered luminosity ranges: for decreasing luminosity, g^m varies from 0.5 to 0.6 
in both the ACS and the WFPC2 samples. Hence, the derived values of ^min are neither 
strictly comparable to one other, nor to the estimates presented in different works. We have 
therefore computed the fraction of binaries with mass ratios larger than a fixed value q = 0.6 
(^>o.6)- This value has been chosen as a compromise between having enough statistics and 
avoiding contamination from single stars (indeed, the line corresponding to g = 0.6 in the 
CMD always runs to the right-hand side of the MS population boundary). In this case the 
"binary population" is made up of stars that, in the CMD, are located between the line of 
constant q = 0.6 (left boundary; see dotted lines in Fig. [3]) and the right-hand boundary 
defined above. Its ratio with respect to the total number of stars gives the fraction of binaries 
with q > 0.6, which is presented in Table [21 Obviously, the obtained values are smaller than 
the corresponding minimum fractions C,min in Table [H We also note that the same behaviors 
discussed above are still present, thus again suggesting that the trend with magnitude could 
be real. One possible explanation for the trend with magnitude could be that bright range 
systematically samples more massive stars, which are also expected to be more centrally 
segregated. 
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4.2. The global binary fraction 



In order to estimate the overall binary content of MIO, independently of the value of g, 
we also computed the global binary fraction (^tot)- This requires us to perform simulations 
of single and binary star populations assuming different input values of the global binary 
fraction (,^in) and then determining ^tot from the comparison between the artificial and the 
observed CMDs: the value of that provides the best match between the two CMDs is 
adopted as the global binary fraction ^tot (see Bellazzini et al. 2002 and S07 for a detailed 
description of the procedure). 

Once we assumed an input value of the binary fraction (^in), for each of the considered 
radial and magnitude bins we have built a sample of A^ms and Ntm stars, with Nbin = A^^in, 
being the number of observed objects (after having taken into accout the number of 
contaminating field stars, discussed in Sect. [3]) in that bin, and A^ms being A^(l — ^mm)- 
The MS stars have been simulated by randomly extracting A^ms values of the mass from 
the present-day cluster m ass function d e rived by Bic|f|, and transforming the masses into 
luminosities by using the iBaraffe et al.l (jl997| ) isochrones. Then, from the artificial-star 
catalogue previously described, we have randomly selected an object with similar (A/ < 0.1) 
magnitude and, if recovered, we assigned its output / and V magnitudes to the considered MS 
star. In order to simulate the binary systems we ran domly extracted N^in values of the mass 
of the primary component from the iKroupa I (120021 ) initial mas s function, arid Nf,in values 



of the binary mass ratio from the f{q) distribution observed by iFisher et al. I (|2005[ ) in the 



solar neighborhood, thus also obtaining the mass of the secondary. After transforming masses 
into luminosities and summing up the fiuxes of the two components, an object with similar 
magnitude was randomly extracted from the artificial-star catalogue and, if recovered by the 
photometric analysis, the shifts between its input and output magnitudes were assigned to 
the considered binary system. Finally, the field stars were added to the sample. The result 
of this procedure is a list of synthetic stars with the same characteristics of real stars and 
containing a given fraction of binaries (^in)- To be precise, the MSs of the resulting artificial 
CMDs are narrower than the observed ones, because the formal photometric errors of the 
artificial star catalogue systematically underestimate the true observational uncertainties. 
This is apparent in Figure HI where, for the magnitude range 19 < / < 19.5, the histogram 



'"'BIO suggested that for stars below O.SM©, the slope of the mass function decreases from 0.23 to —0.83 
moving from the inner to the outer regions (for reference, the slope of the canonical Salpeter mass function 
would be —2.35). In order to understand how the assumed mass function may affect the binary fraction 
estimates, we have re-computed Ctot by adopting the core mass function for the whole cluster. Within 
the errors, the resulting values of £,tot turn out to be in agreement with those presented in Table [3l thus 
guaranteeing that the global binary fraction is just mildly sensitive to changes in the mass function. 
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corresponds to the distribution of the observed color differences A{V — I) with respect to the 
MSRL, the sohd hne is the best-fitting Gaussian of the blue-side of this distribution (gray 
histogram; the red-side has been ignored because it also includes the contribution of binaries 
and blends), and the dashed line is a Gaussian with a dispersion obtained by adopting the 
formal photometric error of the artificial star catalogue. In order to correct for this bias 
and adopt realistic values of the photometric uncertainty, we increased the formal errors ct/ 
and (Ty thus to reproduce the observed error distribution as a function of magnitude. As a 
check, we verified that the width of the resulting color distribution with respect to the MSRL 
well matches the observed one. An example of the synthetic CMD thus obtained, compared 
to the observed one is shown in Fig. [51 From the simulated catalogue we then computed 
the ratio Vsim = ^Mn^/^MS between the number of synthetic stars belonging to the "binary 
population" defined in Sect. 14. 1^ and that of the synthetic "MS population". The same was 
done for the observed data-sets, thus obtaining robs = ^b^/^MS- 

For every value of ^in, from 0.5% to 25% with steps of 0.5%, the entire procedure was 
repeated 100 times. Then, the penalty function x^(6n) was computed as the summation of 
{.fsim,i — fohs)"^ for « = 1, 100, and the associated probability -P(^in) was derived. To illustrate 
this. Figure E] shows the distribution of P as a function of the adopted values of ^in, ranging 
from 3% to 10%. The mean of the best-fitting Gaussian gives the global binary fraction 
(^tot) and its dispersion has been adopted as the error. The values of ^tot obtained in 
the various radial and magnitude ranges are reported in Table [31 The global binary fraction 
shows the same radial behavior observed for ^mmi varying from ~ 14% or ~ 10% in the 
cluster core (for the full and the intermediate magnitude ranges, respectively), down to 
~ 1.5% in the outskirts (for both). As before we find a dependence of the binary fraction on 
the magnitude. This could be an effect of mass segregation, since the average binary mass 
in the bright, intermediate and faint ranges is M ~ 1.1, 0.8, O.5M0, respectively. However 
it could also depend on the assumed mass-ratio distribution and the estimate of blended 
sources, and future studies will be required to resolve this. 



5. DISCUSSION 

We have presented a homogeneous analysis of the binary fraction in MIO as a function 
of the radial distance from the cluster centre, from the core region, out to ~ 2rh. Within 
the errors, the derived core binary fraction is consistent with that measured in other GCs, 
which have typical values of ^tot spanning from ~ 10% to ~ 25% (S07; Davis et al. 2008) 
but it is signi f icantl y smaller than that estimated for the faintest clusters in the sample of 



SoUima et al.l (120071 ). which reach also binary fractions S,tot ~ 50%. This is in agreement 
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with t he quoted ant i -correlation between binary fraction and total luminosity (IMilone et al. 



2008 



radius 



SoUimal l2008l : ISoUima et al.l |2010[ ). Also the binary fraction be yond the half-ma ss 
(~ 1%) is consistent with previous estimates in GCs (see Table l lDavis et al.ll2008l ). 



The minimum binary fraction decreases from ~ 6% within Tc, to ~ 1% beyond the half- 
mass radius. An analogous trend was found for the fraction of binaries with q > 0.6 and for 
the global binary fraction (Fig. [7]), the latter varying from ~ 14% to ~ 1.5% from the core 
to beyond the half-mass radius. Such a radial behavior is in agreement with what has been 
previously found in the few other GCs where this kind of investigation h a s been performed 



{ Rubenstein &: Bailynl 119971 : iBellazzini et al.l 120021 : IZhao &: Bailyru l2005t ISommariva et al. 



20091 . and references in Table 1 of Davis et al. 2008). It is also in agreement with the expecta- 
tions of dynamical models, where the effect is essentially due to the mass-segreg ation process. 



which leads to an increase in the number of binaries in the cluster cores (e.g., iHurley et al. 



20071 : ISollimall2008l : iFregeau et al.ll2009l : Ilvanovall2011l). Indeed , the half-mass relaxation time 
of MIO (~ 0.8 Gyr, Harris 1996; see also iGnedin et al.lll999l : iMcLaughlin fc van der Marel 



2005 ) is just a small fraction (~ 4%) of the cluster age {t ~ 13 Gyr; Potter et al. 201ol ) 



so 



it seems safe to conclude that the system has already had time to achieve equipartition. 



By comparing the radial variation of the MS stellar mass function derived from the 
observations, with that obtained in N-body simulations, BIO suggested that either an IMBH 
or a population of binaries should be present and act as a central energy source in MIO, 
supressing the mass- segregation profile. In particular, the shallow mass-segregation profile 
could be modeled without an IMBH only when the simulations started with a primordial 
binary fraction of about 3 — 5%. Within this framework, in Figure [7] we compare our derived 
values of C,tot, with those obtained from the dynamical evolution of the 5% primordial binary 
population in the 32K particle simulation of BIO. For a proper comparison we considered a 
simulation snapshot at ~ 7 relaxation times, and only those binaries made of two MS stars 
and with the primary component in the mass range 0.44 O.56M0, corresponding to the 
lower and upper cuts of the intermediate magnitude range along the MSRL. The resulting 
binary fractions for the three considered radial bins are: ^N-body = (0.070 ± 0.02), (0.032 ± 
0.007), (0.026 ± 0.006), from the centre to the outskirts. It is apparent from Fig. [7] that 
the observed binary fraction is larger than the simulated one, especially in the core. This 
indicates that the binary content of MIO is indeed sufficient to account for the observed mass 
segregation suppression, with no need to invoke an IMBH as additional energy source. 
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Radial 
bin 


full 
mag range 


bright 
(18.8 < / < 20.3) 


intermediate 
(20.3 < / < 21.5) 


faint 
(21.5 < / < 23) 


r < Tc 
rc<r <rh 
r > Th 


(6.3 ± 0.4)% 
(3.6 ±0.2)% 
(1.2 ±0.3)% 


(7.6 ±0.5)% 
(3.9 ±0.2)% 


(4.6 ±0.5)% 
(3.1 ±0.2)% 
(1.5 ±0.6)% 


(1.1 ±0.3)% 



Table 1: Minimum binary fraction (^mm) of MIO in the three considered radial bins and the 
magnitude ranges defined in a Sect. O 



Radial 
bin 


full 
mag range 


bright 
(18.8 < / < 20.3) 


intermediate 
(20.3 < / < 21.5) 


faint 
(21.5 < / < 23) 


r < Tc 
rc<r <rh 
r > rh 


(5.2 ±0.3)% 
(3.0 ±0.2)% 
(0.8 ±0.2)% 


(6.2 ±0.5)% 
(3.2 ±0.2)% 


(4.2 ±0.5)% 
(2.9 ±0.2)% 
(0.7 ±0.4)% 


(1.0 ±0.3)% 



Table 2: As in Tabled but for the fraction of binaries with mass ratio q > 0.6 (^g>o.6)- 
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Radial 


full 


bright 


intermediate 


faint 


bin 


range 


(18.8 < 7 < 20.3) 


(20.3 < 7 < 21.5) 


(21.5 < 7 < 23) 


r < Tc 


(13.8 ±1.4)% 


(15.1 ± 1.9)% 


(10.0 ±1.6)% 




Tc <r <rh 


(7.4 ± 0.6)% 


(7.6 ± 0.8)% 


(6.3 ± 0.6)% 




r > Th 


(1.5 ±0.6)% 




(1.5 ±1.0)% 


(1.5 ±0.7)% 



Table 3: Global binary fraction {^tot) of MIO for the considered radial and magnitude 
intervals. 
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RA— RAq [arcsec] 

Fig. 1. — Map of the ACS and WFPC2 data-sets. The two circles mark the position of the 
core and half- mass radii {vc = 48", Vh = 147"). The ACS data-set samples the inner portion 
of MIO out to ~ r/i, while the WFPC2 data-set (consisting of the data acquired with the 
three wide field cameras) covers a region between one and two half-mass radii. 
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(V-I) (V-I) 

Fig. 2. — CMDs of the ACS and WFPC2 data-sets. The stars used to estimate the binary 
fraction are plotted in black. The three horizontal lines mark the values along the MSRL 
that separate the magnitude ranges considered in the work (see Sect. E]). The photometric 
errors at different magnitude levels are shown. 
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(V-I) (V-I) 

Fig. 3. — Selection boxes used to define the "MS population" and the "binary population" 
for the two data-sets. The thick solid line marks the MSRL. The dashed lines bound the 
"MS population", made of stars with a color difference from the MSRL smaller than three 
times the typical photometric error at that magnitude level (see Sect. [3]). The gray region 
marks the "binary population" selection box, with its left-hand side corresponding to the 
redder boundary of the MS population region (right-hand dashed line) and its right-hand 
side corresponding to the equal-mass binary boundary shifted to the red by three times the 
photometric error (Sect. 14. ip . The dotted line represents the locus defined by binary systems 
with mass-ratio q = 0.6. 
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A(V-I) 

Fig. 4. — Observed color distribution of MS stars with respect to the MSRL, in the magnitude 
range 19 < / < 19.5 for the ACS data-set (histogram). The sohd hne corresponds to the 
Gaussian that best-fits the blue-side of the observed distribution (gray histogram) , while the 
red-side has not be taken into account since it also includes the contribution of binaries. The 
dashed line is a Gaussian with a dispersion equal to the formal photometric error derived 
from the artificial star simulations (BIO). 
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(V-I) (V-I) 

Fig. 5. — CMDs corresponding to the simulated single MS stars, binaries, and field stars for 
the case of the ACS sample, Tc < r <h, and = 6.3 (upper panels). The lower panels show 
the comparison between the combined synthetic CMD (left-hand side) and the observed one 
(right-hand side). 
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Fig. 6. — Probability distribution of the adopted input binary fractions ^in, for the case 
of the ACS sample, Vc < r < Vh and the intermediate magnitude range. The mean and 
the dispersion of the best-fitting Gaussian give the global binary fraction and its error: 
^r^OT = (6.3 ± 0.6)% (see Table EJ. 
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Fig. 7. — Upper panel: Radial behavior of the global binary fractions (black dots) and of 
the fraction of binaries with mass ratio q > 0.6 (empty dots) estimated for the intermediate 
magnitude range (20.3 < I < 21.5). Lower panel: comparison between the observed values 
of ^TOT (the same as above; black dots), and the corresponding current binary fractions 
obtained from an N-body simulation that started with a 5% primordial value and no central 
IMBH (gray regions; see text and BIO). 



